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Abstract 

The Hultman numbers enumerate permutations whose cycle graph 
has a given number of alternating cycles (they are relevant to the 
Bafna-Pevzner approach to genome comparison and genome rearrange- 
ments). We give two new interpretations of the Hultman numbers: in 
terms of polygon gluings and as integrals over the space of complex 
matrices, and derive some properties of their generating functions. 

Introduction. In the paper [3j on genome comparison and genome rear- 
rangements, Bafna and Pevzner raised the problem of decomposing a permu- 
tation into the minimal number of "transpositions" (here a transposition is 
understood as an exchange of two contiguous intervals of the permutation) . 
An important tool they introduced to deal with this problem is the cycle 
graph of a permutation. We recall that the cycle graph of a permutation 
TT £ Sn, denoted by G(7r), is the directed edge-colored graph with vertices 
{0,1,..., n} and edges of two colors: grey edges going from i — 1 to i and 
black edges going from vr, to 7rj_i, i = 0, . . . ,n (throughout this note we 
assume that vro = and consider i modulo n + 1). An alternating cycle in 
G(7r) is a directed cycle with edges of alternate colors. Notice that at every 
vertex of G^n) there is one incoming edge and one outgoing edge of each 
color. This means that there is a unique disjoint decomposition of the edge 
set of G{tt) into alternating cycles, see Fig. [TJ 

In his thesis [7j, Hultman attempted to characterize the number H{n, k) 
of permutations in Sn whose cycle graph has exactly k alternating cycles. 
These numbers, now carrying his name (see www.oeis.org/Al64652), have 
later been studied by several authors (cf. ^ and [5| to name just few). As 
it is shown in [4J, the Hultman numbers are closely related to the (unsigned) 
Stirling numbers of the first kind S{n,k) (see www.oeis.org/A008275) that 
count permutations in Sn whose disjoint cycle decomposition consists of k 
cycles: 

f if n — A; is odd, , , 

i/(n. A;) = <^ ' (^i^ 

otherwise. 
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Figure 1: The cycle graph G{Tr) of the permutation tt = (2314), where the grey 
edges are drawn by dashed arrows and the black; edges are drawn by solid arrows. 
There are 3 alternating cycles: 0-1-3-4-1-2-0, 2-3-2 and 4-0-4. 

A closed formula for the Hultman numbers was obtained in [S]. 

In this note we give two new interpretations of the Hultman numbers 
in the spirit of [S]: as numbers of certain polygon gluings and as integrals 
over the space of complex matrices. We also give a recursion relation for the 
Hultman numbers and derive some properties of their generating functions. 

Polygon gluings. Consider a 2n-sided polygon, whose boundary consists 
of n black sides followed by n grey sides; the black sides are oriented in the 
counterclockwise direction, and the grey sides are oriented in the clockwise 
direction, see Fig. [2l 




Figure 2: A 27i-gon {n — 4) with n black sides (solid arrows) and n grey sides 
(dashed arrows). The pairs of sides that are glued together by tt = (2314) ^I'S 
connected with dotted lines. 

Pairwise gluing of black sides with grey sides (respecting orientation) 
gives an orientable topological surface without boundary of topological genus 
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g > (the genus g depends on the gluing). At the same time, the boundary 
of the polygon turns into an oriented graph with A; > 1 vertices and n 
edges. The numbers g and k are related by the Euler characteristic formula 
2 — 2g = k — n + 1, so that k = n + 1 — 2g. We denote by hg(n) the number 
of genus g such gluings of a 2n-gon. 

Theorem 1. The Hultman numbers H{n, k) and the numbers hg(n) of genus 
g gluings of a 2n-gon described above are related by the fomula 

H{n,n + l-2g) = hg{n). (2) 



Proof. We start with a slightly different interpretation of the cycle graph 
G{'k). Consider two oriented cycles (that is, 2-regular oriented graphs) of 
length n + 1, one colored in grey and the other colored in black. The vertex 
set in both cycles is {0, . . . , n}, but in the grey cycle the vertices follow in the 
clockwise order, and in the black cycle they follow in the counterclockwise 
order. We identify the vertex vTj of the grey cycle with the vertex i of the 
black cycle (we assume ttq = 0). Obviously, the obtained graph coincides 
with the cycle graph G'(7r), see Fig. [TJ 

We label the black sides of the polygon by numbers from 1 to n in the 
counterclockwise order, and the grey sides by numbers from 1 to n in the 
clockwise order, both times starting from the initial vertex 0. Clearly, a 
gluing of a 2n-gon of the type considered above is uniquely described by 
a permutation tt € Sn, where vTj is the number of the grey side identified 
with the ith. black side. Let us cut the polygon along the diagonal (n,0), 
i.e., we add one black edge and one grey edge connecting the vertex n to 
the vertex 0, see Fig. [21 Now we have two n-gons, one with black boundary 
and the other with grey boundary, whose sides are pairwise identified by 
means of the permutation vr (vro = 0). These two boundaries glued together 
give a graph that we denote by r(7r). The construction is quite similar to 
that of the cycle graph G(7r), but instead of gluing vetices we now glue 
edges according to the same rule. The graphs G{-k) and r(7r) are closely 
related to each other: it is straightforward to verify that there is a one-to- 
one correspondence between the alternating cycles in the cycle graph G{tt) 
and the vertices in the polygon gluing graph r(7r). To complete the proof, 
we recall that k = n + 1 — 2g, where k is the number of vertices of r(7r), and 
g is the genus of the glued surface. □ 

Matrix integral. Denote by M{N) = Matc(A^ x A^) the linear space of 
complex N X N matrices; the (complex) dimension of M{N) is A^^. The 
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space M{N) has a natural Gaussian probabilistic measure 

/ 1 ^ 
diiN=\—-7=^\ e~"^'^^^*^ /\ dxij Adxij, (3) 

i,j=l 

where X = {xij}fj^^ S M{N), the star * denotes the Hermitian conjugation 
and Tr is the trace. Note that the space M{N) equipped with the measure 
liN is also called the complex Ginibre ensemble. 

Theorem 2. Put 

[n/2] 

PniN)=^H{n,n + l-2g)N^-''a+\ (4) 

where H{n, k) are the Hultman numbers. Then 

Pn{N)= [ Tr(X"X*")d/i^. (5) 

Jm{n) 

Remark 1. More general matrix integrals over the space M{N) are con- 
sidered in [Ij. 

Remark 2. Below is a list of the several first polynomials pn{N): 
Po{N) = N, 

P2iN) = N^ + N, 

P3{N) = N'^ + 5N^, 

p^{N) = N^ + 15N^ + 8N, 

p^{N) = N^ + 35N^ + 8AN^, 

Pq{N) = iV^ + 70N^ + 469 A^^ + 1807V, 

P7(iV) = + 126N^ + 1869iV^ + SOAm"^ , 

pg{N) = iV^ + llON"^ + 5985iV^ + 26060iV^ + 8064iV, 

pg{N) = iV^° + 330iV^ + 16401iV^ + 1529007^ + 193248iV^ 



Proof. It is a fairly standard exercise in t'Hooft graphic calculus to reduce 
the matrix integral in Eq. ([5]) to a sum over Feynman diagrams (polygon 
gluings), cf. e.g. [8j, [9j. We will briefly explain how it works. By definition 
we have 

N 

i2n = l 



N 

Tr(X"X*") = X] • • • 

n=l 
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and a simple computation shows that 



M{N) 



XijXkid^iN = SikSji, / XijXkid^N = / XijXkidfiN = 0. 

M{N) J M(N) J M{N) 

Applying Wick's formula (cf. [8], |9]), we get 

Xi^i2 ■ ■ ■ Xi^i^^-^Xi-^i^^ . . . Xi^^^i^^j^duN 

where aj = 2n + 1 — ttj (we assume that i2n+i = ^i)- Therefore, 

■^^W -GS„il=l »2„ = 1 

We note that the pairs of indices {ikik+i} correspond to the black edges of 
the polygon on Fig. and the pairs of indices {ia^+iia^} correspond to 
the grey edges, so there is a one-to one correspondence between the pairings 
of indices and polygon gluings. Moreover, it is not hard to see that for a 
given TT £ Sn 

N N 

n=l «2n = l 

where g denotes the genus of the surface glued from the 2n-gon by means 
of vr. This yields 

. [n/2] 

/ Tr(X"X*") dfXN = E hgin)N-^-'^3+^^ 

Jm{N) 

and Eq. ([5]) now follows from Theorem [TJ □ 



Generating functions and recursions. Here we collect some simple 
facts about the recursive relations and generating functions for the Hultman 
numbers that we did not find in the literature. 
Consider the generating functions 

oo CXD ^ 

Fix,N) = E E H{n,n + 1 - 2g)N^-^3+'^ (6) 

g=0 n=2g 

and 

oo 

Hgix)= ^ H{n,n + l-2g)x^. (7) 

n=2g 
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Theorem 3. We have 

(') 

(ii) H{n,n + 1 — 2g) = hg{n) satisfy the recursion 

{n+2)hg{n) = (2n+l)/ig(n-l)-(n-l)/ig(n-2) + n^(n-l)/ig_i(n-2); 

(Hi) the polynomials PniN) defined by Eq. ^ satisfy the recursion 

(n + 2)pn{N) = (2n + l)Npn-i{N) + (n - l){n^ - N^)pn-2{N) 
with po = N, pi = N'^; 

(iv) 

Ho{x) = j^, H,{x) = j^^^, 9>1, 

where Pg{x) = Ylt^g ^g,i^^ ^■^ polynomial with integer coefficients, 

(2s)! 1 J D /I \ (4g-l)!! 

ag,2g = a<,,49-2 = 1, and Pg{l) = ^^+^. 

Remark 3. Several first polynomials Pg{x) are listed below: 

Po{x) =1, 

Pl{x) =x2, 

P^{x) =x'^(8 + 12x + x^), 

p.^{x) =x^(180 + 704x + 528x2 + 72x^ + x^), 

p^{x) =x®(8064 + 56160X + 98124x2 _^ 53792x^ + 8760x^ + 324x^ + x^), 
P^{x) =x^° (604800 + 6356160X + 19083456x2 + 21676144x^ + 

+ 9936360X'' + 1759520x^ + 103040x^ + 1344x^ + x^). 

Remarkably, all polynomials Pg{x) have positive integer coefficients. More- 
over, the integers -Pg(l) are well known (see www, http://oeis.org/A035319 ) 
- they enumerate genus g orientable gluings of a 2g'-gon [6], or the per- 
mutations in S'4g_i whose cycle graph alternating cycles are all of length 2 

m- 

Proof. Part (i) follows from Eq. ([1]) and the fact that 

oo n „ 
n=0 k=0 
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where S{n, k) are the Stirhng numbers of the first kind. Similarly, the recur- 
sion S{n+1, k) = S{n, k—l)+n S{n, k) for the Stirling numbers immediately 
implies (ii). Part (iii) is a direct consequence of (m). The proof of (iv) is by 
induction on g and follows the proof of Theorem 1 in |2]. The cases = 0, 1 
being easy, assume that the statements of part (iv) of the theorem hold for 
g-1, g>2. Put Hg{x) = x^Hg{x), then the recursion (ii) is equivalent to 
the ODE 

(1 - xfH'gix) + (1 - x)Hg{x) = x^H'^Liix) + 2x^H'^_^ix) 
with initial condition Hg(0) = 0. Therefore, we have 

l-x) 

Jo 

The elementary formula 



Hg{x) = {l-x) I ' (i_^)3 '~'^'' dt. (8) 



immediately yields 



+ 2x' 



(9) 



(l-xfj \{l-x)P^ 

_ Q2(a - + • • • + (/3 - a)2(/3 - a + l)a;"+^ 

- (1 - x)/5+3 

Since, by assumption, 

^9-1 W (l_a,)49-3 (l-x)49-3 ' 

applying Eq. (fTOj) we get that 

x^H'^LAx) + 2x^H'^_,{x) Qg{x) 



(10) 



(1-X)3 (1-^)49+3' 

where Qg{x) = Yli=2g+i 1g,i^^ ^ polynomial with integer coefficients, 
qg,2g+i = {2gf{2g - l)ag-i,2g~2 = 2(2^)!, 

9c/,4g = 2ag_i^4g_6 = 2. 

Consider the Laurent expansion 

4g+3 



(11) 



Qg{x) 



(l-x)49+3 {l-xY' 
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then we have 

l-x ^ i(l-xY 

i=2 ^ ' 

where the initial condition Hg{0) = imphes that 

4g+2 



C = - ^ 

Now put 



=2 



43+2 4g+2 

P^(,) =J2"-^{il- - (1 - = ^ pg,x\ (13) 

i=2 1=0 

By construction, we have pg^o = 0, therefore Hg{x) = Pg{x) / {1— x)'^^'^^ since 
they both satisfy the same first order ODE with the same initial condition. 
Moreover, since hg{l) = ... = hg{2g — 1) = 0, we also have pg^i = . . . = 
Pg,2g+i = 0. Inverting we see that 

ag,2g = Pg,2g+2 = qg,2g+l/{2g + 2) = (2c/)!/ {g + 1), 
03,49-2 = P3,43 = 93,43/2 = 1 

as claimed. We also see that Pg{x) = Pg(x)/x'^ = (1 — x)'^^^^Hg{x) must 
have integral coefficients because Hg{x) does. 

To complete the proof it is sufficient to show that 

(4g-l) (4g-3) (2g-l) 
25 + 1 

(note that Po(l) = Pi(l) = 1). We have 



£v, (l-x)P;_i(x) + (%-3)Pg_i(x) Pg^,{x 



9-iV-^ (l-x)4s-2 (l-x)49-2^ 

(1 - x)P;,(x) + (45 - 2)P,,i(x) _ P^^2(x) 



^g-iy-^) (l-x)49-i {i-xyg-i 

, . _ (l-^)P;2(^) + (4g-l)P3,2(^) 
^3-lW - (1 _ ' 



and from Eq. (|lip it then follows that 

Qg{x) = (1 - x)(x^P;2(a:) + 2x3Pg_2(a:)) + (4g - l)x^P3,2(x). 
From here we easily get 

P,,i(l) = (45-3)Pg_i(l), 

P,,2(l) = (45 - 2)P,,i(l) = (45 - 2)(4g - 3)P3_i(l), 
Q,,i(l) = (45 - 1)P5,2(1) = (45 - 1)(45 - 2)(45 - 3)P,_i(l). 



8 



Clearly, <5g,i(l) = ^3^4^+3 in the Laurent expansion (fT2]) . and from Eq. (fT3]l 
we obtain P,(l) = 3^ = (4g-i)(4g-2)(%-3) ^j^.^^^^ ^ 
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